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Chalcogenide materials have received great attention in the last decade owing to 
their application in new memory systems. Recently, phase-change memories have, in 
fact, reached the early stages of production. In spite of the industrial exploitation 
of such materials, the physical processes governing the switching mechanism are still 
debated. In this paper we work out a complete and consistent model for transport 
in amorphous chalcogenide materials based on trap-limited conduction accompanied 
by carrier heating. A previous model is here extended to include position-dependent 
carrier concentration and field, consistently linked by the Poisson equation. The 
results of the new model reproduce the experimental electrical characteristics and 
their dependences on the device length and temperature. Furthermore, the model 
provides a sound physical interpretation of the switching phenomenon and is able to 
give an estimate of the threshold condition in terms of the material parameters, a 
piece of information of great technological interest. 
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I. INTRODUCTION 



With the recent introduction of the 22 nm node in the fabrication process, the semicon- 
ductor industry seems to be very close to its technological limit. According to the Inter- 
national Technology Roadmap for Semiconductors,- it could eventually be possible to scale 
down the actual devices to the next lx generation only by redesigning either the device, or 
the productive process, or both. It will be more and more complicated, if not impossible, 
to continue this trend much further.-^ Alternatively, new materials have to be explored and 
different architectural solutions implemented. 

In the memory technology, non-silicon materials are being widely investigated in order 
to introduce faster, more scalable, and reliable devices. The present frontier of the research 
is represented by metal-oxide resistive RAMs and conductive-bridge RAMs,- while Phase- 
Change Memory (PCM) prototypes have been studied in the last years^ and are now in 
the early production stage.-^ Phase-change materials like chalcogenides are known from the 
early 1960s 6 and show the property of an easy, reversible transition between crystalline and 
amorphous phases, characterized by significant changes in optical reflectivity and electrical 
resistivity. 7 Due to high optical contrast between the two phases, these materials have been 
employed for optical storage since the mid 1990s. 8 In the new century, the strong difference in 
resistivity characterizing the two phases pushed chalcogenide materials as suitable candidates 
for solid-state nonvolatile memories. Furthermore, some chalcogenide glasses also feature 
an ovonic threshold-switching in the amorphous phase, which implies a negative differential 
resistance (NDR) in the current-voltage characteristic before the phase change takes place. 
Even though the first PCM arrays have already been released to the market, the knowledge 
of the physical process governing the ovonic switching mechanism is still a step behind. The 
availability of a theory able to identify the threshold point and predict the device behavior 
under given operating conditions is still sought by scientists and engineers to tailor the 
materials and the device set up. 

The first microscopic interpretation of the switching behavior was due to Adler and 
coworkers,- 1 ^ who supposed the creation of a micrometer-wide low-resistance filament in 
the amorphous matrix, thus reducing the resistance of the device. Later on, the switch- 
ing behavior was also found for sub-micrometer devices,— and this gave rise to alternative 
interpretations. Using the standard macroscopic quantities of the theory of transport in 

2 



semiconductors, like concentrations, velocities, and mobilities, it is possible to explain the 
switching in terms of drift-diffusion with impact ionization, 12 or cooperative detrapping.— 



Alternative interpretations are due to Karpov and coworkers,— 1 ^ and to Ielmini and 
coworkers.— 1 ^ The interpretation provided by Karpov and coworkers preserves the idea 
of conductive filaments as responsible for the switching behavior. From an energetic bal- 
ance, they derive a model based on nucleation and growth of a thin crystalline filament 
that progressively expands in the amorphous matrix until it connects the electrodes, thus 
dramatically reducing the electrical resistance of the device. 



On the other hand, Ielmini and coworkers proposed a thermally-assisted trap-limited 
conduction mechanism, where the switching is ascribed to the increase in the average kinetic 
energy of the carriers through the device, as a result of the balance between the field-induced 
energy gain and the energy relaxation due to the scattering with phonons. A non-uniformity 
of the electric field is also found. This model has the capability of interpreting not only 
the electrical characteristic of a memory cell, but also its dependence on thickness and 
temperature, which suggests a thermally-activated conduction mechanism. Similarly, it has 
been shown through Monte Carlo simulations that a switching behavior can be triggered 
by a space-charge accumulation near the contacts, if a field-enhanced hopping conduction is 
considered. 18 



According to Ref. Il7| a non-uniformity in the electric field across the device is found near 
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and above the switching condition, which must be sustained by a non-negligib 
charge in the region close to the cathode. However, the equations used in Ref. 
include the effect of the variable concentration of the active carriers across the device, which 
is instead approximated as a constant. In the present paper we use the ideas presented 
in Ref. |l7| as a starting point, and we work out a complete and consistent model which 



includes diffusion and achieves self-consistency between charge distribution and field. The 
results of this theoretical development reproduce the electrical characteristics and their de- 
pendences on the geometrical scaling factors and temperature, and provide a sound physical 
interpretation of the results. 



II. THE MODEL 



The presence of defects inherent to an amorphous material implies the existence of a 
number of localized trap states in the band gap. If the conduction were due to a pure 
hopping process, carriers would tunnel among traps until they reach the collecting con- 
tact. In the so called trap-limited conduction regime, instead, carriers undergo continuous 
trapping-detrapping processes by which they absorb and release energy in such a way that 
they overcome the potential barriers and move across the device. An intermediate mech- 
anism (thermally- assisted tunneling), where carriers absorb some energy and tunnel to the 
next trap, is also possible. The three mechanisms coexist, and the final transfer rate must 
encompass all of them. Since pure hopping is mainly effective in the low-temperature range 
(i.e., well below room temperature),— 1 ^ we can neglect it as we are not interested in that 
temperature range. As for thermally-assisted tunneling, it was shown that its dependence 
upon temperature and field is similar to that of the trap-limited conduction process.— For 
these reasons, we develop the model for the trap-limited conduction case. 

Let ut be the trap concentration, assumed to be spatially uniform. Under equilibrium 
conditions, the concentration of electrons is independent of the position and given by the 
Fermi-Dirac statistics; the electron temperature T coincides with the lattice temperature 
To. In the off-equilibrium conditions, we assume translational symmetry along the x and y 
directions, so that the quantities of interest depend only on z, the current direction. The 
carrier concentration n(z) is still described by a Fermi distribution where a quasi-Fermi 
level Ep(z) replaces the Fermi level Ep Q and T(z) replaces T . A flat trap density of states 
T = n T /AE G is assumed inside the band gap AE G = E c (z) — E v (z), where E c (z) and 
Ey(z) denote the band gap edges. The band gap follows the potential profile along the z 
axis. The number of carriers between Et and Et + dEp is given by: 
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The integration of Eq. ([T]) over the band gap yields the carrier concentration at z: 
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The equilibrium value n is obtained once E c (z), E v (z), E F (z) and T(z) are replaced by 



their equilibrium values Eco, Eyo, Epo and T : 
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When Epq is sufficiently far from the band edges Eco and Eyo and close to midgap, 
no ~ ut/2. Under equilibrium conditions, the material is neutral. This situation may be 
obtained, for instance, by assuming donor-like (acceptor-like) traps and a negative (positive) 
compensating charge no. 

The model is in principle the same if one deals with electrons or holes. For this reason, 
in the followings we would rather term the two leads where carriers enter or leave the device 
as the injecting contact, z = 0, and the collecting contact, z = £, respectively. For the sake 
of simplicity, we develop and discuss the model only for the case of electrons. 

Let Az be the average traveled distance between the sites of successive detrapping- 
trapping events. Following Ref. if the detrapping time is much longer than the traveling 
time, it can be taken equal to the transfer time r. For thermally-activated processes r is 
exponentially dependent on the barrier height experienced by the carriers at the detrapping 
event. Thus, we can define two different times r_> and tv that apply to the motion in the 
two directions, i.e., from z to z + Az and from z to z — Az. The local field gives rise to 
opposite effects, as shown below: 

'Eo(z) -E T + AU{z, z±Az 
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where To is a characteristic transfer time for the process at hand, and AU is the shift of 
the barrier height with respect to the equilibrium value due to the local electric field F(z). 
With reference to Fig. [H let XAz be the distance of the maximum of the energy profile along 
the transition path from the trap having the lower z coordinate, with < A < 1. Using a 
first-order approximation, the shifts AU result: 



AU(z, z + Az) = qF(z)XAz, AU{z, z - Az) = -qF{z){\ - A)Az, 
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with q the absolute value of the electron charge. For the sake of simplicity, A is assumed 
independent of z and equal to 1/2. 

At a given position z, the carrier velocities in opposite directions around the energy Et 
are given by: 
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FIG. 1: (Color online) Schematic view of the energy profile for transitions from z to 
z + Az and from z to z — Az. The dashed profile corresponds to the equilibrium condition, 
while the solid curve is obtained under the influence of a local field F(z) < that enhances 
the transitions towards larger z's. Symbols are explained in the text. 



and, using Eq. (JTJ), their average values over the entire distribution are 
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The current density J{z) thus reads: 
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J{z) 



—n[z) 



(v^z)) - (M*)> 



obtain: 
(9) 



which, using Eqs. p]), (EJ) and ([7]), becomes: 
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If J is fixed, as happens in the description of switching materials whose current-voltage 
characteristics are typically S-shaped, the model requires the determination of three un- 
known functions, namely the electric field F(z), the quasi-Fermi level Ep(z) and the carrier 
temperature T{z). Two additional equations to be coupled to Eq. ffTU]) are then required. 

One of them is the Poisson equation, which, due to Eqs. 02]) and ([3]), reads: 
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where £ = £q£ t is the dielectric constant of the material. In writing Eq. lull) we have taken 
into account that the material under equilibrium conditions must be neutral everywhere. 
Since E F0 has been assumed close to midgap and sufficiently far from the band edges, Eq. 
(II ip simplifies into: 
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The third equation comes from the power balance. Depending on how effective electron- 
phonon scattering is in dissipating the power transferred to the carriers by the electric field, 
the average kinetic energy of the carriers (thus their temperature) may or may not stay 
tied to the equilibrium value. Let AEj° T (z) represent the excess energy, i.e., the difference 
between the actual energy of the carrier distribution in z and the energy that the same 
population would have if kept at the equilibrium temperature T : 
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Here Ep(z) is defined at any z by imposing the constraint 
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which ensures the same population for the two distributions. 

In order to write the power balance, one considers the power flowing through two different 
sections of the device at a distance dz: 
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Here is the energy density flux in z; dip = —F(z)dz is the variation of the electrostatic 
potential in the z direction, and the last term of the r.h.s. represents the power exchanged via 
inelastic electron-phonon scattering. The derivative can be expressed in the relaxation-time 



approximation as^- 
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r r being constant a relaxation time. After expanding the l.h.s. of Eq. ( !T3|) to the first order 
in dz and dividing both sides by dz, one gets 
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The energy density flux $(z) can be calculated following the same scheme adopted for 
the current density J(z) in Eq. (|SJ): 
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represent the two average energy fluxes of the carrier distribution in opposite directions at 
a given coordinate z. 

The r.h.s. of Eq. f TT5"j) can be replaced with its first-order approximation in Az/2, this 
leading to 
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The set of Eqs. ( fTUj) . ( fill and ( fl8l) leads to the determination of the unknown functions 
F(z), Ep(z) and T(z) for any given current density J. Eq. (fl8l) involves the second deriva- 
tives of the unknown functions; nevertheless, a numerical analysis has shown that the term 
proportional to the second derivative in Eq. ffl8l) is negligible with respect to the other term, 
but for a narrow region close to the injecting contact at the highest currents. For the latter 
case this term, though effective to some extent, still remains smaller than the other one 
in the l.h.s. of Eq. ([181) . The second-derivative contribution can thus be neglected without 
substantially affecting the physical results. 

However, an analytical closed form for the solution of the set of Eqs. ffTUl) . (fTTj) and 
(fl8l) cannot be obtained because the integrals in Eq. (fl8l) have to be evaluated numerically. 
One can overcome this problem by replacing the Fermi distribution function with a suitable 
approximation x(Et,z) defined in such a way that: i) x(Ef(z)) = 1/2, ii) x(Et,z) shares 
the same asymptotical values of the original Fermi-Dirac distribution, and iii) the following 
symmetry holds 1 — x{Ef — A) = x{Ef + A). The above requirements are satisfied for 
instance by 
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where a good choice for the parameter Q is Q = 3/4, as shown in Fig. |2j 
Let us introduce the following dimensionless functions: 
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which describe the off-equilibrium local field, the shift of the quasi-Fermi level, and the 
electron temperature, respectively. 

After introducing the x{Et, z) function in place of the Fermi distribution, and neglecting 
the second derivative in Eq. (fl8|) . the set given by Eqs. ( ITUj) . (TTTj) and ( IT8|) can be manipulated 
through the straightforward, though lengthy, calculations summarized in the Appendix to 




Energy (arb. units) 



FIG. 2: (Color online) Comparison between the Fermi distribution function (solid dotted 
line) and the approximating function x{Et, z ) for different values of the parameter Q 
(dashed lines). The best approximation is obtained when Q = 3/4. 



yield: 

^- = N*(g,t) (20) 
dz 

J g (J, 9,t)% + W, 9,t)f z = J* (f, 9, t, ^\ (21) 
H 9 (f, 9,t)^- + H t (f, g,t)^ = H* (f, g, t, (22) 



dz dz \ dz t 

From the above, one immediately obtains: 

d/ =Ap dg _ J*H t - J t H* dt _ J g H* - J*H g 

dz ' dz J g H t -J t H g ' dz J g H t -J t H g ' 1 ' 

The definitions of the symbols can be found in the Appendix. 

Using a first-order Runge-Kutta integration scheme, the equations above are numerically 

solved once the values of f(z),g(z) and t(z) are provided at the coordinate z — 0. Two 

boundary conditions can easily be inferred by supposing that the electrons at the injecting 

contact are at equilibrium, namely Ep(0) = Epo and T(0) = To. The definitions of g(z) and 

t(z) allow for the direct conversion of these boundary conditions into g(0) = and £(0) = 1. 
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The value for f(z) at the boundary must be such that the global charge neutrality of the 
system holds true: 

e \f(0) - F(£)] + q f [n(z) - n ) dz = 0. (24) 
L Jo 

A trial-and-error procedure for the initial value /(0) is applied, until the solutions for f(z), 
g(z) and t(z) satisfy Eq. fl2l]) within the desired precision. 



III. RESULTS AND DISCUSSION 

A. Current-voltage characteristics 

We report in Fig. [3] the I(V) characteristic obtained with the model given by Eqs. ([2D]) . 
f l2~Tj) . and f )22|) (solid black curve), where V is calculated as the integral of the field F(z) along 



the device. For comparison the result in Ref. [17| (dash-dotted red curve) and experimental 



data for a GST-225 memory cell available therein are also shown (dots). These data corre- 
spond to a memory cell with a bottom contact electrode with cross-section S = 1000 nm 2 
and length I = 40 nm. The dashed green curve refers to a calculation performed using 



the present model with the same parameters reported in Ref. [17|. It is clearly seen that, 
for the same set of parameters, the new model would lead to an evident increase in the 
conductivity (the dashed curve yields a given current at a lower voltage) and to a reduction 
of the threshold voltage, without substantially affecting the threshold current. These dif- 
ferences are due to the different approximation adopted for the carrier distribution function 



by the two models. While in Ref. [17| the carriers that contribute to transport are only those 
above the Fermi level, in the present model the Fermi distribution function is integrated 
over the entire band gap, this making the conductivity higher. Since the threshold point is 
determined mainly by the current, as will be discussed later on, the increased conductivity 
implies a smaller threshold voltage. 

In order to fit the experimental data, it has then been necessary to calibrate the param- 
eters. The new parameter set is reported in Table [J along with the set used in Ref. |l7| for 



the sake of comparison. The band gap considered here is compatible with literature data 
for the amorphous GST-225. The two time constants r and r r are consistent with those 
suggested by Mott and Davis for amorphous semiconductors.— 

The agreement between the results of the present model with the new set of parameters 
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FIG. 3: (Color online) I(V) characteristic obtained for the present model compared to 



experimental data and to the corresponding curve taken from Ref. [17|. The dashed green 
curve refers to a "hybrid" configuration where the parameters proposed in Ref. LL7J have 
been used in the present model. The arrows show the positions along the I(V) curve of the 

points Pi. . . P4 cited in the text and in Figs. HI 



parameter 


present work 


Ref. 17 


AE G 


0.68 eV 


0.6 eV 




10 20 cm- 3 eV- 1 


10 20 cm- 3 eV- 1 


TO 


1.2 • 10~ 14 s 


1.0 • 10" 14 s 


T r 


0.78 • 10~ 13 s 


1.0 • 10~ 13 s 


Az 


7 • 10~ 7 cm 


7 • 10~ 7 cm 


e r 


15 


15 


T 


298 K 


298 K 



TABLE I: Parameters used for the best fit reported in Fig. [3j The values used in Ref. 

are also reported for comparison. 
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and the experimental data of is quite good. Differences with the results of model of Ref. 
are found in the NDR region, where the present model estimates larger potential drops for 
any given current. We point out, however, that above the switching point the experimental 
data depend on the characteristics of the external circuit. For this reason they should not 
be considered as for the sake of comparison. In fact, under different experimental conditions 
the current may rise nearly vertically at a holding voltage.- 2 Then, the crystallization of the 
material occurs. Results from the present model obtained far above threshold should thus 
be considered only qualitatively. 



B. Microscopic interpretation of the switching 

In order to analyze the microscopic process leading to the threshold switching, it is useful 
to consider the physical quantities of interest along the device. Apart from the amplitude 
of the variations, a common behavior for the dimensionless electric field f(z) and carrier 
temperature t(z), and for the carrier concentration n(z), can be outlined. The calculated 
profiles for these quantities are reported in Fig. H]for the four points Pi. . . P4 shown in Fig. 

m 

Since the injecting contact acts as an infinite reservoir of carriers in thermal equilibrium 
with the lattice, carriers enter the device at thermal equilibrium. However, due to the im- 
balance between the power provided by the field and the power loss due to electron-phonon 
scattering, they tend to heat up as long as they travel along the device. This effect is neg- 
ligibly small at the lowest currents. The increase in the carrier temperature would enhance 
the flux towards the collecting contact and, in order to keep J at its prescribed value, it is 
compensated by a decrease in the carrier concentration and in the electric field, consistently 
with the Poisson equation. Thus, a positive charge accumulates near the injecting contact. 
The position of the quasi-Fermi level shifts towards a lower energy to account for such a 
change in the carrier concentration. 

Going farther from the injecting contact, the heated carrier population enables a more 
effective dissipation through inelastic scattering, which prevails over the power provided by 
the electric field. The electron temperature, after reaching a maximum value, decreases. 
A thermal overshoot is thus created near the injecting contact. Since the charge must 
vary continuously in space, the presence of a positive charge implies that the electric field 
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FIG. 4: (Color online) Dimensionless electric field f(z) (top), dimensionless electron 
temperature t(z) (middle), and electron concentration n(z) (bottom) for the four currents 
identified in Fig. |3j Note the scale change between the OFF and the ON regions at 

z = 7 nm. 
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continues to decrease also in the region where the carrier temperature is reduced, although 
at a lower rate. In order to preserve the current, the carrier concentration must compensate 
the reduction in velocity due to the smaller field and temperature, so that it increases 
significantly and approaches again the equilibrium value. When the charge neutrality is 
attained far enough from the injecting contact, the electric field and the carrier temperature 
saturate as well. It follows that a saturation value for the position of the quasi-Fermi level 
is also found. 
As in Ref. 



16l we refer to the OFF region as the zone of length £off close to the injecting 
contact where the physical quantities vary most appreciably, and to the ON region as the 
longer zone of length Iqn where they have their saturation values, as shown in Fig. HJ 

Even though the interpretation given above applies at any current, in the first part of 
the characteristic (point Pi, I = 5 nA) the electric field is low and does not provide enough 
power to determine an appreciable electron heating. This fact implies that the electron 
temperature and concentration are tied to their equilibrium values across the whole device 
(Ohmic behavior). 

A similar situation applies also for the most part of the subsequent exponential region. 
However, as the current increases, the electric field f(0) grows rapidly (P2, / = 0.5 /j,A) 
and eventually triggers an appreciable electron heating. In turn, the latter causes a slight 
depletion of the OFF region and a corresponding weak non-uniformity of the electric field. 
Due to the exponential relationship linking the current and the field, this picture is more 
and more evident as the switching current is approached (P3, 1 = 4 yuA). 

As the current is increased above the threshold point (P4, I = 30 fiA), /(0) still continues 
to grow, thus inducing a larger carrier heating in the OFF region. As a consequence the 
local electric field rapidly falls, and the electron concentration is strongly depleted to keep 
the current constant. The high value of the electron temperature in the OFF region reflects 
into a high saturation value also in the ON region, which is the key condition to restore the 
equilibrium between energy gain and energy relaxation. Since in the ON region the carrier 
concentration has always the equilibrium value at any current, every further increase in the 
carrier temperature can only be compensated by a further decrease in the electric field. For 
such a reason, the resulting electric field is smaller than that found at the switching point. 
Since the ON region is substantially much longer than the OFF region, a smaller field in the 
ON region leads also to smaller potential drop across the device. 
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Finally, in the region above the switching point the power density dissipated and trans- 
ferred to the lattice via electron-phonon scattering is high. This may give rise to lattice 
heating, which is a favorable condition for the creation of a local crystalline nucleus that can 
eventually evolve into a crystalline filament. The incorporation of the Fourier heat equation 
and a local lattice temperature in the model goes beyond the scope of the present paper and 
is planned for a future work. 

In conclusion, few words must be spent on the role of the Poisson equation, by comparing 
the present results with those of Ref. [l7L As in the ON region charge neutrality is kept, 
both models provide similar results; on the contrary, in the OFF region the introduction of 
a self-consistent non-uniform carrier distribution yields a more accurate physical picture. In 
particular, the presence of a minimum in the carrier concentration profile implies a thermal 
overshoot and an initial different curvature of the field profile. 

C. The switching condition 

The last part of this section is devoted to the analysis of the switching condition as a 
function of the device length (Fig. [5]) and the lattice temperature (Fig. [6]). 

The OFF region, as shown in Fig. HI extends over about the first 5 nm from the injecting 
contact, which is the space where carriers cannot fully relax the power provided by the field 
through electron-phonon scattering and heat up. The microscopic phenomena occurring 
in the OFF region suggest that £off must primarily depend on material properties, like 
the density of traps tit, the position of the Fermi level with respect to the bottom of the 
conduction band, and the relaxation time r r , but must be independent of the device length. 
As a consequence, when the latter exceeds approximately 2£off the potential drop in the 
ON region dominates over that in the OFF region, and the threshold voltage scales almost 
linearly with the device length, as shown in the inset of Fig. |5j In this case, the error made by 
considering the electric field F t h in the ON region as representative of the field in the entire 
device at threshold is negligible. On the other hand, small deviations from linearity are 
found for shorter devices, as the potential drop in the OFF region gains relative importance 
over that in the ON region. According to the present model, ultra-short devices are not 
expected to show a NDR portion of the I(V) characteristic, but, rather, a steep rise of the 
current with an almost constant potential, in agreement with the results of Fig. [5j 
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FIG. 5: (Color online) /(V) characteristics for devices with different lengths ranging from 
10 nm to 100 nm with steps of 10 nm. The inset shows the linear relationship between the 

threshold voltage and the device length. 



In order to assess the effect of the lattice temperature on the switching condition, we 
report in Fig. [6] the I(V) characteristics for different lattice temperatures. It is found that 
the threshold current increases with temperature, while the threshold potential decreases, as 
typical of chalcogenide glasses.— Moreover, the calculations show that the threshold poten- 
tials tend to accumulate towards a limiting value as the temperature increases, suggesting 
the existence of a minimum threshold electric field. 

With the purpose of understanding this phenomenon, we recall first that the ON region is 
defined as the zone where the equilibrium carrier concentration is restored, and the electric 
field and carrier temperature saturate. This physical condition is expressed by making the 
l.h.s. of Eq. ( TlTl) to vanish; viz., 

Taking F t h as representative of the field within the device, the product JthFth is the input 
power density w t h at the threshold point. Here and in the following, the suffix th indicates 
the quantities evaluated at the threshold point. The definitions of A_E™ T (z) and U(t,a,(3) 
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FIG. 6: (Color online) /(V) characteristics for different lattice temperatures in the range 
198 K to 498 K. The intersections with the dashed lines represent the switching points as 

calculated from Eq. (12] 



given by Eqs. (lA.4j) and (1A.9D in the Appendix yield: 

Ith &<fth 



w th 



t 2 th U(t,a,P) -E/oK>$) 

r(kT ) 2 th 



th 



If tth < 2.5, as happens for a wide range of lattice temperatures up to 600 K (see also Fig. 
Ej), the exponents in U(t,a,(3) can be neglected, and the above equation further simplifies 
to 

1 a T 2 — T 2 
w th = ^-Tk 2 ^-^, (25) 



where the condition that the quasi-Fermi level in the ON region is close to the local equilib- 
rium value (gth ~ 0) has also been used. The lack of an analytical solution for t(z) makes 
it impossible to further simplify Eq. (1251) into a compact form showing the dependences on 
the parameters of the material only. However, by means of numerical analyses it is possible 
to obtain an empirical formula linking tth with the lattice temperature. 

Since the threshold current strongly increases with temperature overcompensating the 
decrease in the threshold potential, the input power at threshold increases indefinitely. The 
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above considerations suggest an approximate function for t t h like, e.g., 



Hh 



6T - T* 

Tn — T* 



(26) 



where the parameter T* has the physical meaning of an asymptotical temperature up to 
which the NDR regions of the I(V) characteristics are possible. For the case at hand, by 
best fitting the data reported in Fig. [7] (dashed green line), one finds 8^2/3 and T* « 730 K, 
a value slightly larger than the glass transition temperature of the material. We point out, 
however, that these values should be considered with CcirG, clS they depend non-linearly on 
the parameters of the model. Numerical calculations have shown that a major role is played 
by the energy relaxation time r r , which significantly influences the threshold current and 
voltage. In fact, increasing the energy relaxation time by one order-of- magnitude reduces 
the asymptotical temperature by about 100 K, whereas reducing the relaxation time by one 
order-of-magnitude let the asymptotical temperature raise only by 35 K. The 9 coefficient 
ranges instead from 0.75 down to 0.63 for the same variations of r r . 
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FIG. 7: (Color online) Values at the threshold point of the dimensionless electric field and 
of the carrier temperature. The dashed green line is the best-fit approximation of tth 

obtained by means of Eq. ( )26l) . 
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By inserting Eq. (|26|) into Eq. (125]) . one finds: 

U,th= 9r r A£ G pT^F ( } (27) 

that represents the critical power at threshold as a function of the lattice temperature. In 
order to provide an estimate of the switching point, this formula must be combined with 
Eq. (1A.2j) that expresses the current in the ON region. Taking again into account that the 
quasi-Fermi level in the ON region is close to its local equilibrium value, and using Eq. ( |26|) 
to express the carrier temperature at threshold, it is possible to obtain the threshold current 
as a function uniquely of the dimensionless field f t h- By inserting the expression for J t ^ from 
Eq. ( 1A.2I) into Eq. ( 1271) . after some algebra, one gets: 

f th smh(f th ) = -— exp — — ) kT — — — — - — - (28) 

9 r r \2kT J (kT -kT*) 2 Q th (T ,T*,6) 

where Q t h(To,T*,0) is calculated by setting t = t t h and g = in Eq. ( 1A.8I) . Once f t h is 
known from Eq. ( 128]) . the threshold potential and current are given as Aip t h = F t h£ and 
J t h = Wth/Fthi respectively. The above results allow for a rough estimate of the threshold 
point from the physical parameters of a given switching material. 



IV. CONCLUSIONS 

An enhanced model for transport in amorphous chalcogenides has been worked out to 
achieve self-consistency between the electric field present along the device and the local 
density and energy distribution of carriers. This characteristic was, in fact, missing in the 
models available in the literature. 

The model relies on the hypothesis that conduction can be described by means of trap- 
limited transport, i.e., sequences of detrapping events letting a carrier move above the 
conduction-band mobility gap, followed by fast recaptures by different traps. Within this 
framework, it is possible to write two equations for the charge and the energy fluxes where 
two characteristic times are present: the detrapping time and the energy relaxation time. A 
third equation expresses charge-field consistency. The carrier concentration and the charge 
and energy fluxes are obtained by means of the integration of the Fermi distribution function 
over the band gap, which is the energy region filled by trap states. The solution of the above 
set of equations yields the three main physical quantities that represent the unknowns of the 
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problem, namely, the carrier concentration, the carrier temperature, and the electric field as 
functions of the position along the device. 

After a proper identification of the physical parameters, the calculated current-voltage 
characteristics are in good agreement with experimental data. Moreover, the model correctly 
identifies both the temperature and device-length dependences of the /(V) curves. 

The ovonic switching behavior of chalcogenide glasses is explained by the effect of carrier 
heating. In trap-limited conduction, conductivity is determined mainly by the detrapping 
time, which depends on the energy barrier separating the energy of the carrier sitting in a 
trap from the conduction-band mobility edge. Detrapping processes are therefore favored 
by a lowering of the barrier due to the electric field and by the raising of the carrier energy 
due to their heating. At and above threshold, any further increase in the current requires 
a higher carrier temperature, which reduces the field in the largest part of the device, so 
that the potential drop decreases. The model can thus predict the switching point starting 
from the physical parameters of the material, which is a key achievement for technological 
applications. 
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Appendix: Analytical derivation of the constitutive equations 

This appendix contains in some details the calculations which transform Eqs. (flOl) . (|TT|) 
and (|18p in the main text into Eqs. (1201) . (1211) and (122]) so that the interested reader can 
follow the mathematical development between the two sets of equations. 

Let us introduce the dimensionless auxiliary variables: 

a (z) = C + g(z), a (z) = ^\, (3 (z) = C - g(z), and f}{z) = -r 

t{z) t(z) 

with Co = AEc/2kTo. Following the calculations sketched in Sect. [Til after replacing the 
Fermi-Dirac distribution with x{Et, z \ one finds (omitting the indication of the dependence 
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on z of the unknown functions): 

»E C 



n 



/ T X (E T ,z)dE T = TkT tN(a,(3) 



(A.l) 
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kT n 
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' 2kT 
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TOT 
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7"0 



cosh(/) exp(-/3 )Q(t, a , a, /3 , 0) 



X(E T , z) - x{E T , z) (Et - E v ) dE T 
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T(kT ) 2 [t 2 U(t,a,f3)-U (a ,f3' ) 

Ec Az 
Tx(E T ,z)(E T - E v ) — exp 
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2TAz{kT n ) 2 
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kT n 
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sinh(/) exp(-P Q ) S(t, a , a, A), /3) 



£c Az 
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' 2kT 



T() 



cosh(/) exp(-p )S(t, a , a, A), /3) 
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iV(a, /3) = a + — exp I — -a — - exp I — -/5 



Q(t,a ,a,Po,P) 



9 - 16t 2 



exp(-a ) 
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At 
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(A.5) 
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(A.6) 
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4t 
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9 43^2 _ g 

S(t, a , a, A,, /?) = %_ 16t 2 + 9 ^T^ip + ex P(~ a o) 



1 / At 



2 V At + 3 



exp | --a 



4i 



4t-3 



2f 



At -3 



exp ( A) - ~p ) . 



(A.10) 



The U (a' , (3' Q ) function in Eq. (1A.4j) can be calculated by setting (t, a, (3) = (1, a' Q , f3' ) in 
U(t, a, (3), with a' and (3' fulfilling the constraints tN(a, (3) = N(af , (3' ) and a' + (3' = 2C . 
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Under equilibrium it is / = 0, g = 0, and t = 1. As a consequence, a eq = (3 eq = Co- The 
equilibrium electron concentration thus reads: 



n = TkT C 



By means of this relationship and of Eq. (lA.lj) . Eq. ffTTj) can then be recast in terms of the 
dimensionless auxiliary variables as follows: 



Tz~ Cl 



tN(a, P) - C 



(A.ll) 



with d = q 2 TAz/(4e). 

Eq. ( flOl) can be approximated by means of Eqs. ([9]), ( 1A.2I) and ( 1A.3I) . Recalling that 



da dp da 1 / dp dt 

— a 



dz dz dz t V dz dz 
and observing also that 



d/3 dp 



dz dz ' dz 



dp 1 /dp dt 



t V dz ^ dz 



dQ 9 , , 4dg 
— = 2exp(-a ) + 
aa 3 aa 



and 



<9g _49Q 



Eq. (jTD|) can be recast as 
dp dt 



Jfl dz + J< dz = sinh (-^(*' a °' a ' ^ 



A_d/ 
Az dz 



+ C 2 Jexp(/3 



(A.12) 



where 



J 3 = cosh(/) 



Q(t, ao, a, /3 , P) + 2 exp(-a ) + - + 



4 i\ <9g /4 i\ ag 



+ 



3 t J da \3 t J dp 



J t = cosh(/) 



'9Q _ a^g _ 0dQ* 
~~dt ~ l~da~ ~ l~dp 

and C 2 = 2T /[qTkTo(Az) 2 }. 

Similarly, after neglecting the second derivative as indicated in the text, Eq. (I18j) can be 
recast making use of Eqs. f lT7|) and ( 1A.5[) . this yielding: 



9 dz dz 



exp(/3 ) 2C 2 Jf - C 3 t 2 U(t, a, P) - U {a', P') 



H f f z , (A.13) 



with 



H g = sinh(/) 



Hf = cosh(f)S(t, ao, a, Po, P), 

S(t,a ,a,Po,P) + 7 I « aa J + 



t \da dp J da dp 
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H t = smh(/) 

where C 3 = 7"o/(4AzT r ). 

Eqs. f lATTj) . flAJ2|) and (1A7L3]) can be turned into Eqs. (120]) , (j2T|) and (f22l by lett 



AT* = Ci 



tN(a,P) - C 



2 d/" 

J* = sinh(/)Q(t, a , a, /3 Q , /3) ( — - — ) + C 2 J exp(/3 ), 



and 



iT = exp(/3 ) |2C 2 J/ - C 3 [t 2 U(t, a, /3) - U (a', /?')] } - H 



dz' 
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